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Abstract 

We analyze the production of p-mesons in p + A reactions including both the 
production by proton-nucleon as well as pion-nucleon collisions within a coupled 
channel transport approach. The final state interactions of the p-meson with 
nucleons are evaluated from a resonance model which allows to extract elastic 
and inelastic cross sections. We include the latter final state interactions, the p- 
meson decay into two pions as well as the final pion-nucleon interactions within 
the transport approach. We find the invariant mass distribution of pion pairs to 
be sensitive to the p-meson properties in the nuclear medium. However, due to 
the strong final state interactions of pions only light targets like ^^C might be 
suited to extract a p-signal from the uncorrelated two pion background which 
carries information about the in-medium properties of the p-meson. 



* Supported by Forschungszentrum Jiilich 
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1 Introduction 



The restoration of chiral symmetry at high baryon density or temperature nowadays is 
a question of primary interest. As suggested a long time ago the short hved p-meson 
is a primary candidate to study possible in-medium mass shifts or dispersion relations 
as e.g. predicted by Brown - Rho scaling or QCD sum-rules 0. Whereas the 
dileptonic decay of the p-meson appears as an ideal probe of the in-medium properties 
of the p-meson, the two pion hadronic decay in principle carries the same information 
provided that the decay products suffer only little rescattering. 

At SPS energies the dilepton data of the CERES and HELIOS-3 Collaborations 
[^] have indicated a possible mass shift of the p-meson following the analysis of Li, Ko 
and Brown [Q as well as Cassing et al. f^, ^ |^. A recent summary of the situation 
has been given by Drees . However, the approaches in and ||^ differ in the actual 
assumptions on the p-meson mass shift. In the model of Ref. |^ a dropping of the 
p-mass at normal nuclear matter density po of about 30-35% is employed whereas in 
the transport study of Ref. [|, ^ || the p-mass drops only by about 18% at po in 
line with the QCD sum-rule study of Hatsuda and Lee Such drastic in- medium 



effects should also be seen in proton-nucleus or pion-nucleus collisions |T0|, |Tl[ via the 
electromagnetic decay of the p-meson or even the a;-meson in case of tt + A reactions 

In the present paper we explore the possibility if such medium effects of the p-meson 
might also be measured by two-pion invariant mass distributions in proton-nucleus re- 
actions. Our study thus is organized as follows: In Section 2 we present the elementary 
p-meson production cross sections from nucleon-nucleon and pion-nucleon collisions 
and determine the elastic and inelastic cross sections of the p-mesons with nucleons 
from a resonance model. In Section 3 we discuss possible in-medium modifications 
of the p-meson as incorporated in the present simulations. In Section 4 we test our 
transport approach with respect to pion spectra in proton-nucleus reactions. Section 
5 is devoted to the actual p and pion dynamics in p + ^^C reactions at 2 - 2.5 GeV as 
well as to the reconstruction of the p-mass distribution from two-pion invariant mass 
spectra, while Section 6 concludes our work with a summary and discussion of open 
problems. 



2 Elementary cross sections 
2.1 p-meson production 

In a previous study we have investigated the production of vector mesons from proton- 
nucleus reactions from subthreshold energies to a few GeV employing empirical spectral 
functions |T^ . It was found that at subthreshold energies the dominant contribution to 
p-meson production stems from secondary vr -|- interactions while at energies above 2 
GeV the primary proton- nucleon channel becomes more important. In the latter study 
the cross section for p-meson production from vr + collisions was separated into two 
parts: i) an exclusive cross section related to binary processes [irN — >■ pA^), which is 
dominant at energies close to the reaction threshold, and ii) the inclusive production of 
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p-mesons at higher energies which is dominated by phase space and can be described 
by the LUND string model [|l3l. 

The cross section for the exclusive reaction ir + N ^ p + N can be obtained - in the 
limit of vanishing interference terms - via an incoherent summation over all baryonic 



resonances decaying into pN channels M, [T5|, |I6| as 



AT J- pN 



where Mr and Tji are the mass and total width of the baryonic resonance, while F^^r 
and Fpjv are the partial width for the initial and final state, respectively. Here k is 
the pion momentum in the resonance frame of reference. We account for a momentum 



dependence of the resonance width as [|Tj, |T3] 



(2) 



where kpt is the momentum at i/i = Mr, / is the minimal orbital angular momentum 
of the system and 6"^ = 0.1 GeV/c. 

The properties of the baryonic resonances are available for resonances up to a mass 
of about 2 GeV PDI. In order to describe the exclusive tt + iV — ^ p + N reaction at 
energies above ^/s ~ 2 GeV we incorporate additionally two 'effective' resonances with 



masses of 2.4 and 3.3 GeV; all resonance properties adopted are listed in Table.|l|. 



Table 1: Properties of baryonic resonances used in the calculations. Apart from the 
known resonances we have added the 'effective' resonances El and E2 to describe 
properly the exclusive production of p mesons in nN reactions. 



Resonance 


Ir 


Mr (MeV) 


Tr (MeV) 


-^Nn (%) 


^Np (%) 




5/2 


1680 


125 


60 


5 




3/2 


1700 


100 


10 


10 


Pu 


1/2 


1710 


110 


15 


20 


Pl3 


3/2 


1720 


150 


10 


70 




7/2 


1990 


260 


20 


5 


D33 


3/2 


1700 


280 


15 


30 




1/2 


1900 


150 


10 


45 


F35 


5/2 


1905 


350 


10 


72 


P31 


1/2 


1910 


250 


25 


37 


El 


3/2 


2400 


1000 


70 


25 


E2 


5/2 


3300 


1000 


60 


8 



+ p ^ p^ + p 



181 are shown 



The experimental data for the exclusive reaction vr^ 
by full squares in Fig. |l|a) as a function of ^/s — ^so, with i/i denoting the invariant 
energy of the colliding particles while the threshold energy is given by y/so = rriN + rnp. 
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Figure 1: Cross sections for p-meson production in vr"'" +p (a) and p + p interactions 
(b). The solid lines show the results from the LUND string model for the inclusive 
production cross section at higher energies, while the dashed line in (a) is calculated 
within the effective resonance model (|l]). The dotted line in (a) is the parametrization 

and show the cross sections for 



from [0. The experimental data are taken from 
the exclusive reaction vr+p — > p+p (full squares) (a) and pp — > p'^pp (full squares) (b) 
as well as for the inclusive reactions vr+p — >■ p° + X (full circles) (a) and — p° + X 
(full circles) (b). 



The calculated cross section within the resonance model (|I|) is shown in Fig. [I|a) by 
the dashed line. Note that for energies above — ^/sq ^ 0.3 GeV a reasonable 
reproduction of the experimental data is only possible due to an incorporation of the 
two 'effective' resonances El and E2 as given in Table 1. The dotted line in Fig. [l|a) 



shows the parametrization from Ref. |]12[ for comparison and practically coincides with 
our present result from 

Since the available experimental data show no difference |ll9| between the reactions 
vr~ + p ^ + n, 7r~ + p ^ P~ + P and vr^ + p ^ p~*" + p we assume in the following 
that the cross sections are the same for vr"*', tt" and vr^, for neutrons and protons as 
well as for p"*", p~ and p°-mesons, respectively. 

At low energies above threshold only the exclusive reactions discussed above are 
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energetically allowed. However, an additional production of pions, i.e. vr + — 
p + A^ + tt's, is possible above ^/so > m^. Since there are not enough experimental 
data available for the inclusive p-meson production from ir + N collisions that allow 
to construct reliable parametrizations, we calculate the cross sections for the inclusive 
reactions within the Lund-String- Model (LSM) from [|I^ . The calculated results within 
the LSM can be fitted by 

ain + p^ p + X) = 3.Q[mb] x (x - l)^'^^ x x'^'^^, (3) 

where x = s/sq and Sq = {rriN + mp)^. Our fit to the LSM results is shown by the solid 
line in Fig. |l]a) and reasonably well reproduces the available inclusive data on p-meson 
production (full circles) at higher energies. 

For the p-meson production from nucleon-nucleon collisions we use the parametriza- 



tion obtained within the LSM as 1112 



a{p + p ^ p + X) = 2.2 [mb] x {x-iy-^' x x''-^ (4) 

with X = s/sq and Sq = {2mN + ^p)^ for the nucleon-nucleon channels. The fit to 
the LSM results is shown by the solid line in Fig. |l|b) in comparison to the experimental 
data for the inclusive production of p- mesons (full circles) from For comparison 
we additionally show the cross section for the exclusive reaction p + p — > p + p + p° (full 
triangles) from |jl8|, which about 1 GeV above threshold is already small compared to 



the inclusive cross section. Furthermore, we assume that the cross sections are isospin 
independent for all possible nucleon-nucleon channels allowed by charge conservation. 

Fig.|I] clearly demonstrates that the cross section for p-meson production from pion 
induced reactions is significantly larger than that from nucleon-nucleon collisions at the 
same available energy above the threshold. As a consequence the pion induced channels 
dominate for p meson production also in p -|- A reactions close to the absolute threshold 
energies as shown in Ref. |l! 



2.2 pN interactions 

The absorption of p-mesons in the nucleus due to the reactions p + N imr + 
with m > 1 substantially attenuates the cross section and spectrum of p-mesons in 
p + A collisions. To account for the p-meson interaction in finite nuclei we need the 
p + N cross section at nuclear densities from to po ~ 0.17 fm~^. Since detailed 
Brueckner calculations for these reactions in the medium are not available we assume 
that the in-medium cross sections are approximately the same as those in the vacuum. 
In the latter case the total and elastic pA^ cross sections can be calculated within the 
resonance model as 

where M^, and Fp^v denote the mass, total and partial width {R pN) of the 
baryonic resonances, respectively. In Eq. (5) k now is the momentum of the p-meson 
in the resonance frame. The summation is performed again over all the baryonic 
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Figure 2: Cross sections for p + p interactions: a) shows the total cross section; 
the sohd hne is the calculation with the resonance model (j^), while the dashed line 
indicates the quark model relation The solid circles show the experimental data 
extracted from vector dominance (p!TD; the squares are the data from [^|. b) The solid 
line shows the elastic cross section from the resonance model while the dashed line 
results from the cr-exchange calculations (|^). The dotted line and the triangles indicate 
the cross section for the reaction p + — vr + obtained by detailed balance (OjO). 



resonances listed in Table.|T| thus neglecting the interference between the resonances. 
The p+N total cross section is calculated with Tout = while the elastic cross section 
is evaluated with Tout = ^pN- 

The result of our calculations for the p + N total cross section is shown in Fig.^a) 
(solid line), while the elastic cross section is shown in Fig.^) (solid line); both are 
displayed up to y/s — rriN — nip ~ 0.3 GeV, only, because we do not include in our cal- 
culation the heavier baryonic resonances El and E2 due to possibly improper branching 
ratios. In view of the neglect of interference terms and the limited energy range con- 
sidered our calculations should only be considered as an estimate that will have to be 
improved in future. We note, however, that a momentum dependent width according 
to Eq. (|^) compared to a fixed average width does not change these cross sections in 
Fig. I sizeably. 
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Following the quark model relations one can express the cross section for the p + N 
interaction at high energies in terms of pion-nucleon cross sections as 



(6) 



The total p+N cross section according to Eq.(|^) is shown in Fig.^) by the dashed line. 
For the vrA^ cross sections entering here we have used a Regge fit to the experimental 
data from |3^. Note, that there is a reasonable agreement between the results from 
the resonance model and the quark model in the overlapping energy region. 

The calculation of the p + N p + N elastic scattering in an effective meson- 
exchange model may be restricted to a diagram with the cx-exchange ||21|. In the 
t-channel the elastic cross section then can be expressed as 



1 m 



a 



pN^pN 



N 



Stt 4s 



dx M{s, x) 



(7) 



where s is the square of the total center-of-mass energy, is the nucleon mass and 



the isospin averaged squared amplitude is given by [|2T], P2| 



M{s,x) 
where 



mi 



1 - 




2 24 ) <^^N '^pp 



{tun + nipY] [s - {niN - m^)^] 
2s 



X - 1] 



In Eq.(R|) we use the monopole form factors 



aNN 



A? 



^^2 



Al + q' 



crpp 



(9) 



(10) 



The cut-off parameters Ai = 0.9 GeV and A2 = 1.0 GeV as well as the constant 
glpp/4:7r = 14.8 were taken from the recent analysis of Friman and Soyeur p3 
meson photoproduction off nucleons, while the coupling constant g^j^^/A-n = 
rrirT = 550 MeV were taken from the Bonn potential [^]. The elastic cross section 



on p- 
8 and 



calculated within the a-exchange model is shown in Fig.^D by the dashed line and 
reasonably agrees with that from the resonance model (solid line). Though the actual 
magnitude of this cross might change by about 30% in a more sophisticated analysis, 
it becomes clear that p — N rescattering will be very important for the final p-meson 
spectra. 

Furthermore, in a vector dominance model the photoproduction of a p-meson is 
related to the p + N cross section as 



2 "yp 647r dcT^p—^pp 



dt 



fill 



where a is the fine-structure constant and the coupling constant 7p/4'7r = 0.61 accord- 



ing to Ref . . The solid circles in Fig.|^a) show the total p + p cross section calcu- 
lated from the experimental forward p-photoproduction amplitudes from PH EBI, E^. 
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Furthermore, the squares in Fig.|^a) show the p + N cross section extracted from p- 
photoproduction on the deuteron (7 + (i — > p° + (i) independently of the vector domi- 
nance model 



The cross section for the p+N n+N reaction channel may be obtained exploiting 
detailed balance as 

1 ql 

and is shown by triangles in Fig.pD) using the experimental cross section of the inverse 
reaction from Ref. [0]. 

Note that Eq. (|l^) is only valid in case of stable particles. For a broad resonance as 
the p-meson it is easy to show p9| that 



where g/ is given as 



2 _ [s - (m^ + wy] [s - jm^ - wy] 



and rUp, Tp are the mass and width of the p-meson, respectively. The cross section for 
the reaction cXp+p^Tr+p calculated by Eq.(|l^) within the prescription (|1]) for cr,r+p^p+p is 
shown in Fig.0b by the dotted line and differs only slightly from the prescription (p!2D . 

Note, that the results obtained from quite different approaches are in reasonable 
agreement with each other. At low energies Tp < 200 MeV the contribution from the 
inelastic reaction channels to the total cross section are dominant and much larger than 
the elastic cross section. Starting from a p kinetic energy of about 300 MeV the total 
cross section then stays constant at roughly 30 mb. Following the quark model the 
inelastic reaction channels to the total cross section are dominant here, too. 



3 In-medium modifications of the p-meson 

In the literature there are several approaches on the in-medium modifications of the p- 
meson. Without going into a detailed discussion we simply adopt two different models 
to explore the possibility if such modifications can be observed in the two-pion decay 
channel. 

Following the results of Hatsuda ||30[ obtained from QCD sum-rules the in-medium 



mass m* of the p-meson can be approximated as 

m* = niB {1 - apN/po) (15) 

where a = [0.16 ±0.06], 7/1^=0.7699 GeV is the bare mass and pn is the actual nuclear 
density, while po = 0.17 fm~^. 

Recently, the p-meson mass in finite nuclei was calculated by Saito, Tsushima and 
Thomas within the quark-meson coupling model. The ratio m* /itlb from ^ 



for C is shown in Fig.p by the solid line as a function of the radial distance r. The 
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Figure 3: The ratio of the in- medium to the bare p-meson mass in a ^^C-target 
as a function of the distance r from the center of the nucleus. The sohd hne shows 
the numerical results from the dashed line is our result employing the linear 

parametrization (|15]) in a local density approximation for a parameter a = 0.16 in Eq. 
(|15|) while the dotted line is obtained for a = 0.14. 
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dashed line in Fig.|^ illustrates the density dependence of the p-meson mass according 
to (^) for a = 0.16 using the ^^C density distribution from |3l|, The dotted line is 
obtained for a = 0.14, respectively. Note that the local density approximation works 
quite well in comparison to the calculations from Ref. pTI, |3^. Indeed, the actual 



numerical results in |^T], ^ show only a weak deviation from the linear parametriza- 
tion [0, 33, 34 1 for a fixed nuclear density. In view of the uncertainty in the parameter 



a it thus appears justified to adopt the local density approximation furtheron. Accord- 
ingly, in the following calculations we will use the scaling (^) while keeping the width 
of the in-medium p-meson equal to the free one. 

Furthermore, according to the hadronic model calculations of Klingl and Weise ^5, 
or Rapp, Chanfray and Wambach |37, 3^ or Friman and Pirner |39] both the mass 



36 



and the width of the p-meson are modified in nuclear matter. In our present study we 
use the results from ^5|, Actually, for the p — >■ tttt decay mode we need the general 
spectral function An = dMn/dM^,.^, i. e. the distribution of the resonance mass Mr 
with respect to the invariant mass of the two pions M^^r- The latter is related to the 
quantity R^^^ from as 



A 



R 



9p 



m 



B , 



677^ 



X R 



1=1 



(16) 



where gp = 6.05 is the strong coupling constant, rriB the free p-meson mass and ^/s the 
invariant mass of the two pions. The solid line in Fig.^ shows the spectral function ([T6| ) 
of the p-meson at density pn = Po calculated with R^^^ from Ref. . 

For numerical purposes we have parameterized (0) at invariant energies above 2 

in the form 

Ar = , r- \o . « (17) 



with the parameters 



2ro 

VllB (1 - 0.34pAr/po) 
0.5 + 0.5pAr/po, 



Fq = 151 MeV and being the bare mass of the p-meson. The fit (|1^) is shown in 
Fig.^ by the full dots and reasonably reproduces the spectral function from Ref. |^5|, ^ 
at densities po/2 < pa? < po. 

The dashed line in Fig.^ represents the p-meson spectral function in free space as 



A 



R 



27r(v/i 



+ r2/4 



(19) 



with rrip = tub while Tp is the momentum dependent p-meson width (^, respectively. 
Furthermore, the dotted line in Fig.^ shows the spectral function (|1^) with the free 
p-meson width and an average in-medium mass rUp = m* from Eq. (|l^) for a=0.16. 

In the following calculations we will use the different spectral functions in Fig.^ for 
a study of the p-meson production in p + A collisions accounting for the in-medium 
modifications as discussed above. 
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Figure 4: The spectral function of the p-meson from Khngl and Weise ^ at 
normal nuclear matter density po (solid line) in comparison with the spectral function 
in free space (dashed line); the dotted line displays the p-meson mass distribution 
according to Hatsuda and Lee at po for a = 0.16 employing the free width of the 
p meson. The full dots show our parametrization ([T7|) for the spectral function from 
Refs. 0. 



11 



4 Background from uncorrelated pions 



A proper calculation of the pion spectra from p + A collisions is especially important 
for the study of p-meson production via the two pion decay mode (p nn), because 
multiple pion production from the reactions p + A imx + X with m > 2, which 
are not associated with p-meson decays, substantially contribute to the invariant mass 
spectrum of two pions and have to be subtracted. 

We use the coupled channel transport model [O] for the description of p + A 



reactions that allows to account for the final state interactions of all hadrons and 
especially for those of the pions from the p ^ tttt decay in the medium. The transport 
approach |^ is known to describe experimental data on pion production in a wide 



range of bombarding energies Here, we briefly present a comparison with the 

experimental data available in the energy region of interest. Since the semiclassical 
transport approach involves a local density approximation explicitly, this also 
serves as a test for the underlying assumptions of the dynamical model itself. 

Fig.^ shows the spectrum of Tr^-mesons produced in p + ^^C collisions at 2.1 GeV 
at the pion emission angle of 2.5° in the laboratory. The data (full dots) from Ref. 



4^ are shown as a function of the Feynman variable xp = pI /Pmax^ where is the 
longitudinal momentum of the pion in the center-of-mass system of the incident proton 
and nucleon at rest while p*max is defined as 

Pmax = ^^7^ • (20) 



The solid line shows the scaling behaviour of the cross section from which describes 
p + A reactions from a bombarding energy of 1 GeV up to very high energies. The 
histogram represents our calculations within the transport approach and shows a good 
agreement with the experimental results. 

Fig.p shows the invariant cross section for 7r"^-meson production from p + ^^C col- 
lisions at 2.5 GeV at an emission angle of 40° in the laboratory. The experimental 



data from the KaoS Collaboration [^] are shown by the full dots while the transport 
calculations are presented in terms of the histogram. We slightly overestimate the pion 
spectrum at high momenta and underestimate for p^^ ^ 0.5 GeV/c. These deviations 
are basically due to the local density approximation employed in the transport model 
which is more uncertain in the surface of a light nucleus such as ^^C. Though keeping 
in mind these minor deviations we conclude that the pion production and the proton- 
nucleus dynamics are described quite well within the transport approach such that we 
now can proceed with the p-meson and two-pion dynamics. 

Due to the different kinematics for the production of multipion states and p-mesons 
from p + A collisions the distribution of pion pairs in momentum space should be much 
broader; thus kinematical cuts could be used in order to enrich the contribution from 
p-meson decays. In this respect we show in Fig.^ the distribution in longitudinal and 
transverse momentum of pion pairs from p + ^^C collisions at 2.5 GeV, that do not 
stem from p-decays; Fig. displays the corresponding momentum distribution of 
pion pairs from p-decays, respectively. As expected, the phase space of pion pairs from 
p-meson decay is substantially smaller than that for the uncorrelated pions. Using the 
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Figure 5: Invariant cross section for vr "'"-meson production from p + ^^C collisions at 
2.1 GeV at an emission angle of 2.5° in the laboratory as a function of the Feynman 
variable xp- The full dots show the experimental data from [O]; the solid line is the 



scaling from while the histogram shows our result calculated within the transport 
approach. 
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Figure 6: The invariant cross section for vr "''-production from p + C collisions at 2.5 



GeV at the pion emission angle of 40° in the laboratory. The experimental data [45 



are shown by full dots while the histogram represents the result from the transport 
calculations. 
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Figure 7: The longitudinal and transverse momentum distribution of pion pairs from 
p + ^^C collisions at 2.5 GeV; a) shows the results calculated for uncorrelated pion pairs 
while b) is for the pion pairs from p-decay. The rectangles indicate the cuts used in 
our analysis. 



kinematical cuts as shown in Fig.^ {p^ < 2 GeV/c, pt < 0.8 GeV/c) one can slightly 
enhance the p-meson signal in comparison to the background. In the following analysis 
we thus will always adopt the kinematical cuts < 2 GeV/c and pr < 0.8 GeV/c. 



5 p-meson production from p + A collisions. 

The inclusive production of p-mesons for p + ^^C reactions has been calculated in 
Ref. [T^ on the basis of empirical spectral functions for energies from the absolute 
threshold up to a few GeV. Since in this work we consider bombarding energies above 
2 GeV the semiclassical transport approach employed should yield similar inclusive 
production cross sections when employing the same cross sections for the pN and vrA^ 
production channels, respectively, and discarding any final state interactions of the 
p-mesons. In fact, by comparing the inclusive p" production cross sections in this limit 
from our transport approach with those from |]12| for bombarding energies of 2 - 3 



GeV (cf. Table 0) we find that both models differ only by about 20 %. The latter 
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deviations between the two approaches are due to the local density approximation and 
the on-shell assumption in the transport model. Since these differences are lower than 
the uncertainty in the elementary production cross sections we can use the coupled 
channel transport approach quite confidently for the following analysis, however, can 
include additionally all final state interactions of the p-mesons as well as for the pions 
which are more important in view of their large cross sections with nucleons. 



Table 2: Comparison of the inclusive p°-production cross section in mb for p + 
reactions within the spectral function approach (SFA) [|l^ and our coupled channel 
transport approach (CBUU). For this comparison all final state interactions of the 
p-mesons have been discarded. 



T (GeV) 


SFA 


CBUU 


2.0 


0.36 


0.43 


2.5 


1.2 


1.5 


3.0 


2.6 


2.9 



5.1 p-meson dynamics 

We first investigate the 'clean' p'^-meson production from p + A collisions by neglecting 
the background from uncorrelated two pion pairs. The spectra of tt^tt^ pairs are 
shown in Fig.^ for p + ^^C and p + "^^^Ph at a bombarding energy of 2.5 GeV as a 
function of their invariant mass. The calculations are performed for a bare p-meson 
spectral function (Fig.|^) and with an in-medium modification of the p-meson (Fig.^) 
according to Eq.dTSD for a = 0.16. Note that the spectra for p + ^"^C are displayed by 
hatched histograms and are scaled by a factor of 10 in order to compare the results for 
the light and the heavy target more directly. 

The low energy part of the invariant mass distribution reflects the pion pairs from 
p-decays inside the nucleus that have suffered strong final state interactions. Obvi- 
ously, such events are dominating for the heavy nucleus ^^^Pb. Even without medium 
modifications the p-resonance is practically melted and the mass distribution is close 
to a thermalized spectrum. For ^^C the rescattering of pions from p-decays is signif- 
icantly reduced and a 'bare' p-meson signal or an in-medium p-meson signal can still 
be identified. In Fig.^) we additionally show the calculated results for the invariant 
mass distribution of the pions from p decay for a ^^C target employing the parameters 
a = 0.16 and 0.14, respectively, as in Fig. 3. The direct comparison demonstrates 
the relative sensitivity to the parameter a in the dropping mass scheme of Eq. (15); 
furthermore, it also demonstrates the relative uncertainty arizing from the local density 
approximation in the transport approach since in view of Fig. 3 this uncertainty can 
be simulated by a corresponding change of the a parameter. 

This initial study shows that the pion mode has substantial disadvantages com- 
pared to the dileptonic p-meson decay not only for the investigation of the in-medium 
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Figure 8: The invariant mass distribution of pion pairs from p-decays for p + 
(hatched histogram) and p + "^^^Pb colhsions at 2.5 GeV calculated with the bare (a) 
and in-medium (b) spectral function of the p-meson according to Eq. (|15|) for a = 
0.16; c) shows the calculated result in the dropping mass scenario for a = 0.16 (solid 
histogram) and a = 0.14 (dashed histogram), respectively. The spectra for are 
scaled by a factor of 10 in a) and b). 
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Figure 9: The distribution of the p-decay rate as a function of the nuclear density for 
p +^'^C colhsions at bombarding energies of 1.5, 2.0 and 2.5 GeV. 

properties of the p-meson but also for a conventional analysis oi p + A ^ p + X re- 
actions that rely on the reconstruction of the p-meson via two pion invariant mass 
distributions. The lightest target nuclei, however, may be used if the invariant mass 
spectrum is measured with a sufficiently high accuracy. 

Of fundamental interest for the p-meson dynamics - relevant for in-medium mod- 
ifications - is the distribution of the decay rate as a function of the nuclear density. 
Obviously, this distribution depends on the beam energy as well as on the size of the 
target and should be properly calculated within the transport approach. In this re- 
spect we display in Fig.^ the number of p-mesons decaying at the nuclear density pn 
for p + ^"^C collisions at bombarding energies of 1.5, 2.0 and 2.5 GeV. All distributions 
are normalized to 1 in order to illustrate the relative variation with the beam energy. 

Note that events with densities < 0.05 fm~^ can be associated with a p-meson 
decay at the periphery of the target or in the vacuum and carry no information about 
the medium effects. With increasing bombarding energy the nucleus becomes more 
transparent for the propagation of p-mesons due to lower total pN cross sections and 
due to the fact that more p-mesons decay in the vacuum. Thus low energy collisions 
are more suitable to study the medium modification of p-meson properties, however. 
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Figure 10: The p- meson decay rate as a function of the nuclear density for different 
targets at a bombarding energy of 2.5 GeV. 



one should keep in mind that the cross section for p-meson production at low energies 
is substantially smaller. 

Fig.|10|, furthermore, shows the p-decay rate as a function of the nucleon density at 
2.5 GeV for the targets ^^C, ^"^Al, ^^Cu and ^ospft. While for the P6-target a major 
fraction of p-mesons decays at normal nuclear matter density, this fraction decreases 
with the target mass significantly. However, even for ^^C still 50% of the p's decay at 
densities 0.5po-po such that modifications of the p-spectral function in the medium can 
still be tested with ^"^C. 



5.2 Calculations with a bare p spectral function 

In actual experiments the invariant mass distribution of pion pairs from p-decay as 
presented in Fig. 8 is swamped by the background from pion pairs that stem from the 
reactions pN — >■ AA — >■ 2ttNN and pN — >■ RN 2'kNN , where R is an intermediate 
resonance that decays into the two pion channel (cf. Fig. 7a). The actual invariant mass 
distribution of pion pairs to be expected experimentally (within the cuts < 2 GeV/c 
and pt < 0.8 GeV/c) is shown in Fig.|ll]a) for p + at 2.5 GeV when neglecting 
any medium modifications of the p-meson. The solid histogram in Fig.|TT]a represents 
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Figure 11: The invariant mass distribution of two pions from p + collisions at 
2.5 GeV using the bare spectral function of the p- meson for < 2 GeV/c and pt < 
0.8 GeV/c. In a) the solid histogram shows the total vr — vr invariant mass spectrum 
while the dashed histogram indicates the p signal after subtraction of the combinato- 
rial background. The solid line is a Breit-Wigner distribution with the bare p-meson 
properties normalized to the total p-production cross section from the transport ap- 
proach as expected in case of no final state interactions. In b) the histogram shows 
the p-meson signal from a) while the dotted line is the statistical spectrum; the dashed 
line is a Breit-Wigner function with the p-meson properties Mq = 766.9 MeV and 
To = 173 MeV and the solid line shows the sum as discussed in the text. 



the pion pair mass distribution including the background pions as well as those from 
p-decay while the solid curve indicates the contribution from p-decay as obtained from 
the total p production cross section without any final state interactions. The latter 
solid line thus would represent a maximum p signal if all p-mesons would propagate 
out of the nucleus and decay in vacuum which, of cause, is unrealistic. 

In order to subtract the background from uncorrelated pion pairs we adopt a combi- 
natorial method, i.e. mix pions from different events. This procedure was demonstrated 
in a recent experimental study of A(1232)-resonance excitation in heavy-ion collisions 
by Badala et al. to be very effective for the extraction of a true correlation signal. 
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The dashed histogram in Fig. 11a) shows the difference between the real and mixed 
invariant mass spectra and corresponds to pion pairs from p-decays including all final 
state interactions of the p-meson as well as its decay pions (cf. Fig. 8a). As discussed 
above the deviation of the solid line in Fig. 11a) from the dashed histogram is only due 
to p and pion final state interactions. Consequently, the p spectral function appears 
substantially distorted by conventional hadronic interactions. 

According to Figs.|^,^ a sizeable fraction of the p-mesons decays at very low nuclear 
density for ^^C or in the vacuum where its spectral function is well known. We thus 
can describe this 'outside' component approximately by a Breit-Wigner function 

(21) 



dM^^ (M,, - Mo)2 + ro74 

where Cq is some constant which (in absence of p-meson melting) is related to the total 
cross section of p-production from p + collisions while Mq and Fq are the mass and 
width of the bare p-meson, respectively. 

The TT-TT mass distribution, that is associated to the p-meson production from p + 
^^C collisions (after background subtraction), is shown again in Fig.[Tl|b) by the solid 
histogram. Apart from the 'outside' decay contribution a further fraction of p's decays 
'inside' the medium at nonzero density. Tentatively, we analyze this 'inside' component 
(together with all effects from final state interactions) in terms of a 'thermodynamical 
approach'. In this limit one can describe this part within the statistical model similar 
to the formalism developed in ||4^, ^ by 



-iVe.pf-[^--^°n, (22) 



dM^^ V ^0 

where the normalization A^, average excitation energy Eq and dispersion Tq can be 
fitted to the p-mass spectrum for lower invariant mass M^^^^ as shown in Fig.|lT|b) by 
the dotted line. Accordingly, the solid line in Fig.p^p) shows our fit to the p-meson 
signal by the sum of the Breit-Wigner function and statistical spectrum (p^), where 
the first one involves the parameters Mo=766.9 MeV and Fo=173 MeV and is shown 
by the dashed line in Fig.pn]b). 

So far, we have still discussed the results calculated with the bare p-meson proper- 
ties. The final state interactions of p-mesons and pions from the p — ^> tttt decay yield 
an invariant mass distribution which is close to the recent predictions from Klingl and 



Weise |^ as discussed in Section 4 (cf. also Fig. 4). Any additional in- medium effects 



thus should show up in a further distortion of the invariant mass distribution. 



5.3 In-medium modifications 

As discussed above the QCD sum-rule approach by Hatsuda and Lee [|] predicts a shift 
of the p-mass with density according to Eq. (^) which corresponds to a shift of the 
p-meson pole. In order to explore if such an effect could be seen in the two pion decay 
mode we have performed calculations for p + ^"^C at 2.5 GeV employing Eq. (|15|). The 
solid histogram in Fig.[12|a) shows the total distribution in the invariant mass of two 
pions while the dotted histogram indicates the background obtained by mixing of pions 
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Figure 12: The invariant mass distribution of two pions from p + at 2.5 GeV 
employing the in-medium modification of the p-meson according to Hatsuda and Lee 
[§] for a = 0.16. In a) the sohd histogram shows the total tt—tt invariant mass spectrum; 
the dotted histogram displays the background from uncorrelated pion pairs while the 
dashed histogram gives the p-mass distribution after subtraction of the background. 
In b) the solid histogram shows again the p-meson signal from a); the dashed line is a 
fit with Eq.(p2D while the hatched histogram is the difference between the p signal and 
the fit Eg. (|2^) , which can be described again by a Breit-Wigner function (solid line). 



from different events. The latter is practically the same as that for those calculations 
performed with the bare p-meson properties (cf. section 5.2). The dashed histogram in 
Fig.[T^) represents the p-signal as obtained by background subtraction. It is seen that 
not only the shape of the mass spectra, but also the absolute normalization are different 
from the calculations with the bare p-properties. The p-meson signal is shown again 
in Fig.|l^) in terms of the solid histogram, which at first glance is quite complicated 
to analyze in order to extract in-medium properties of the p-meson. 

We proceed as in section 5.2 and fit the spectrum with the function (E^) as shown 



in Fig.|T^) by the dashed line. The hatched histogram in Fig.|T2|b) then shows the 
difference between the p-mass spectrum and distribution (P^); it can again be fitted by 



a Breit-Wigner distribution (|2T|) with M =708 MeV and r=270 MeV. In comparison 
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with Fig.|TT]b) the dropping of the p-meson mass thus appears possible to be extracted 
from the invariant mass distribution of two pions at least for light nuclei as ^^C. 

We have also studied the modification of the p spectral function following the model 



from Klingl and Weise For simplicity of the analysis we have 'switched-off ' 

the final state interactions of pions stemming from /9-decay in order to distinguish the 
change of the spectral function due to such processes from the in-medium modifica- 
tion (15). Fig.|TB| shows our results calculated for p + ^^C collisions at 2.5 GeV. The 
hatched histogram stems from p-mesons, which decay outside or at the periphery of 
the nuclear target, while the solid histogram illustrates the pion pairs from the 'inside' 
decay of the p-mesons. Note that now the internal component has no pronounced res- 
onance structure anymore and is similar to the spectrum from pion rescattering. The 
dashed histogram in Fig.^ shows the sum of these two components and the solid line 
is a fit (0) with the parameters M =765 MeV and r=240 MeV. The 'extracted' mass 
of the p-meson in this case is close to the bare mass while the width is considerably 
larger. 



6 Summary 

Within the coupled channel transport approach we have studied the production of p- 
mesons in proton-nucleus collisions and the detection of p-mesons via the pion decay 
mode (p — > TT -|- tt) employing different spectral functions for the in-medium p-mesons. 

The elementary cross sections for p-production from p + N and vr -|- interactions 
as well as total and elastic p + N final state interactions were calculated within the 
resonance model at low invariant energies and in quark-string models at higher invariant 
energies, which provided consistent results in a reasonable agreement with available 
experimental data. Note, however, that we have discarded interference effects in the 
resonance model and that more elaborate approaches will be necessary in future to pin 
down these cross sections more accurately. 

Special attention was paid to the background from the production of two uncorre- 
lated pions as well as for the attenuation of the p-meson yield due to pN interactions 
as well as the final state interactions (vr -|- A^) of pions emerging from p-decay. We 
found that only light nuclear targets might be suitable for the experimental study of 
p-production via the pion decay mode. Appropriate kinematical cuts and background 
subtractions in principle allow to reconstruct the mass and width of the p-meson. How- 
ever, for such a study a good mass resolution for the pion 4-momenta is required as 
well large detector acceptance. 

We have investigated different schemes for the modification of the p-meson in nu- 
clear matter and its consequences for the invariant mass distribution of two pions that 
can be detected asymptotically. The vr — vr invariant mass spectra calculated within 



the hadronic scenario of Klingl and Weise [^, ^ are substantially different from those 
obtained by simulations based on the QCD sum-rule approach by Hatsuda and Lee 0. 
These differences are much larger than the theoretical uncertainties arizing from the 
local density approximation in the transport approach and might actually be studied 
experimentally by large area pion spectrometers e.g. at COSY, CELSIUS or the SIS. 
However, according to the authors opinion the actual shape of the p-meson spectral 
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Figure 13: The invariant mass distribution of pion pairs from p + at 2.5 GeV 
including the in-medium modification of the p-meson according to Khngl and Weise 



50]. The hatched histogram shows the spectrum from p-mesons decaying 'outside' the 
nucleus, the solid histogram is the 'inside' component, while the dashed histogram is 
the sum of both. The solid line shows a Breit-Wigner fit with M= 765 MeV and r=240 
MeV. 
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function cannot be determined accurately enough from the data analysis. 
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